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Minghao Wang/’^ El L- F- Wei/’^ and J. Q. Liang' 

‘institute of Theoretical Physics, Shanxi University, Taiyuan 030006, China 
^Quantum Optoelectronics Laboratory, School of Physics Science and Technology, 

Southwest Jiaotong University, Chengdu 610031, China 
^State Key Laboratory of Optoelectronic Materials and Technologies, 

School of Physics and Engineering, Sun Yat-Sen University, Guangzhou 510275, China 

(Dated: February 16, 2015) 

The Berry phase (BP) in a quantized light field demonstrated more than a decade ago ( Phys. Rev. Lett. 89, 

220404) has attracted considerable attentions, since it plays an important role in the cavity quantum electrody¬ 
namics. Flowever, it is argued in a recent paper ( Phys. Rev. Lett. 108, 033601) that such a BP is just due to 
the rotating wave approximation (RWA) and the relevant BP should vanish beyond this approximation. Based 
on a consistent analysis we conclude in this letter that the BP in a generic Rabi model actually exists, no matter 
whether the RWA is applied. The existence of BP is also generalized to a three-level atom in the quantized 
cavity field. 

PACS number(s): 42.50.Pq, 03.65.Vf, 42.50.Ct. 


I. INTRODUCTION 


Thirty years ago Berry discovered that, 'when a quantum 
system varies slo'wly around a closed loop in a certain parame¬ 
ter space the eigenstate of Hamiltonian 'will acquire a geomet¬ 
rical phase in addition to the usual dynamic phase 101. Shortly 
after, this phase, called usuaUy Berry phase (BP), was gener¬ 
alized to various versions |lll[3>[lall31- Now, it is found that 
geometric phases are related to many physical problems such 
as Aharonov-Bohm effect, quantum Hall effect, and Born- 
Oppenheimer approximation, etc. H. Hopefully, geometric 
phases play an important role in fault-tolerant quantum com¬ 
puting Jl5i|l8|] . 

The most BPs inyestigated preyiously are in the semiclassi- 
cal context, namely, the quantum systems are driyen by a clas¬ 
sical field. In 2002, Fuentes-Guridi et al. J?!] generalized the 
original Berry model of a spin-1/2 particle in a time yarying 
magnetic field to a full quantum counterpart, wherein the clas¬ 
sical driying field was replaced by a quantized field. A photon- 
dependent BP is then found in the usual Jaynes-Cummings 
(JC) model. Interestingly, such a BP still exists eyen when the 
field is at yacuum. The existence of the BP ifoll in the two- 
leyel atom system has been successfully extended to yarious 
models, including multi-atom Dicke model [0, [I^ |2^ and a 
multileyel atom in a quantized field ifTH] . 

Surprisingly in 2012 Larson claimed that, the BP in quan¬ 
tum optical system is just a result of the rotating waye ap¬ 
proximation (RWA) and should yanish beyond this approxi¬ 
mation 10]. Note that this argument is obtained based on a 
semiclassical Rabi model, wherein the field operator is simply 
replaced by a complex C-number, and thus the releyant inter¬ 
pretation is practically not related to the quantized light field. 
Alternatiyely we present, in this letter, a uniyersal formulation 
of the BPs for the Rabi model and show that the non-zero BP 
always exists, no matter whether the RWA is applied. 


‘Electronic address: sesiseu@gmail.com 


The general description of the BP in a Rabi model is giyen 
in Sec.II. Then the yalidity of such a BP is yerified by inyesti- 
gating the semiclassical counterpart of Rabi model in Sec. III. 
In Sec. IV we generalize our results to a three-leyel atom in 
the quantized light field beyond the RWA. Finally we present 
the conclusion and discussion in Sec.V. 


II. GENERAL EORMULATION OE BP IN RABI MODEL 

The interaction of a two-leyel atom with a single-mode 
quantized field can be generally described by the Hamiltonian 


Hq = uja^a + —CTz + X{a+a + a-o') + + (T_a), 

( 1 ) 

in which and a± = + iay)/2 are called the pseudo¬ 

spin operators for the two-leyel atom of the eigenfrequency v. 

(a) denotes bosonic creation (annihilation) operator of the 
single-mode quantized field with the frequency w. Obyiously 
when Ajvr = 0 ,namely under the usual RWA, the Hamilto¬ 
nian reduces to that of the JC model. Correspondingly, when 
A NR = A it becomes the standard Rabi-model Hamiltonian. 

Following Fuentes-Guridi et al. |0] the BP can be obtained 
in terms of an unitary transformation |0] 

U{g)) = exp^—iipa' a) (2) 

applied to the Hamiltonian Hq such that 


H{(p) 


U{ip)HQU\ip) 

uja ^a' + —CTz + \{a+a' + a-a 
+XNR{cr+a'^ + a^a'), (3) 


with a'{(p) = a exp (/</:). The eigenstates |'0n,(fp)) of the 
Hamiltonian H{<f) are obtained as \tl}n{yp)) = (7(i^) ]■!/)„) 
with \tj)n) being the eigenstates of Hq. When the angle yari- 
able tp slowly yaries from 0 to 27r a BP giyen by 
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7„ = i j) (tpni'f)] ^ \'>Pn{‘f)) = 27r a+a |^/;„) (4) 

is generated for the eigenstates |r/)„((^)). When Xnr = 0, Hq 
reduces to the JC model Hamiltonian with the ground state 
denoted by |r/)o) = |0) 0 \g)- Here, \n) is the Fock state of n 
photons for the field, and \g) (|e)) is the atomic ground (ex¬ 
cited) state. It is seen from Eq. (|4| that the BP is zero for the 
ground state {ipo), but non-zero for any excited state even if the 
filed is at the vacuum. One can easily prove that the non-zero 
BP is practically the half of the solid angle subtended by the 
traversed loop of the eigenstate !'!/)„((/?)) on the relevant Bloch 
sphere in the basis of |n -f 1)0 Ig) and |n)0 |e). Here, the half 
solid angle is 7r(l — cos0„), with being the angle between 
the eigenstate vector and the north axis. Note that is also 
associated with the coefficient of the eigenstate, and thus the 
above solid angle can be expressed as 2Tr[{ipn \ |'0„) — n]. 

On the other hand, when Xnr = X i.e. the case of Rabi 
model, any eigenstate can no longer be written as the form 
of |0) 0 {Cg \g) -f Ce |e)) due to the existence of the counter 
rotating wave terms. This implies that, the average photon 
number of any eigenstate of the Rabi model should not be 
zero. As a consequence, the BP of any eigenstate induced by 
the unitary transformation U{(p) is always non-zero according 
to the Eq. (|4]l. In fact the non-zero BPs have been found for 
various eigenstates in the Rabi model ifl^ . 


III. APPARENT CONTROVERSY ON BP IN T HE 
SEMICLASSICAL THEORY OE RABI MODEL AND ITS 
RESOLUTION 


In this section we first of all briefly analyze how the claim 
of vanishing BP comes out for the Rabi model, and then 
show that it is actually non-zero if the proper semiclassical 
approximation is made. 


A. Vanishing BPs as a result of improper semiclassical 
approximation 

Following Ref. d, we begin with the ’’semiclassical ap¬ 
proximation” of the Hamiltonian Eq. Q (quotation mark here 
means improper) i.e., 

Hcip) = uj\a\ + -cr^ 

+X{ae^^a+ + a*e-^‘^a_) 

-I-ae*‘^(T_), (5) 

which is obtained by directly replacing the bosonic operators 
a and with the complex C-numbers a, a* respectively. It 
can be further written as 

Hcip) = uj\af+ \a\ cos (j){X +XNR)crx 

+ \a\ sin (j){XNR - X)ay+ ^az, (6) 


with a = |a| exp(i(/?') and = p + ip'. Obvi¬ 

ously this Hamiltonian is formally equivalent to that of 
a spin-1/2 particle driven by a magnetic field; B = 
(|a| cos0(A -I- Aath), |a| sin(/)(AAr_R — A), v/^). The 

eigenstates of the effective Hamiltonian can be found as a 


spin coherent states where ft ■ a 


±nj = ± 



with ft being the unit vector along the i3-direction. There¬ 
fore the slow variation of the parameter p from 0 to 27r cor¬ 


responds to the eigenstates 


± n ) traversing a loop on the 


Bloch sphere and then the eigenstates will acquire BPs given 
by 7± = ±H/2, where H is the enclosed solid angle by the 


loop. On the other hand if the eigenstates 


± n ) traverse just 


an arc but not a closed loop, the relevant geometric phase is 
zero Sin. 

Specifically when Xnr = 0 the Hamiltonian Hc{p) be¬ 
comes the ’’semiclassical approximation” of the JC model. Its 
eigenstates are found as 


\L+) 


/ cos f \ 
Vsinfe-*'?^/ 


|i-) 


f 

Vcos 


( 7 ) 


where tan6* = 2 |a| A/A, A = v — oj. It is seen that, 
when p varies slowly from 0 to 27r the eigenstates \L±) 
acquire the BPs 7 ± = ±H/2 with the solid angle H = 
27r(l — cosO). While for the ’’semiclassical approximation” 
of the Rabi model when A nr = A, the ground state becomes 


\R-) 


1 /2|a|cos0A 

^E_{2E_ - v) \E- - v/ 2 ) ’ 


( 8 ) 


with E- = w \a\^ — /A + 4 \a\^ cos^ </>. Since the ele¬ 
ments 2 |a| cos 4> and E- — j^/2 ^ 0 are both real, the ground 
state |i?_) traverses only an arc giving the zero BP when p 
varies from 0 to 27r. Thus the conclusion of vanishing BP in 
the Rabi model d is recovered and it seems that the non-zero 
BP in the JC model is just a result of the RWA. 

However, this apparent controversy comes entirely from the 
improper semiclassical approximation used in Ref. d- The 
Hamiltonian Hc{p) by directly replacing the creation (anni¬ 
hilation) operator (a) with a complex C-numbers a* (a) 
does not correspond to the original one H{p) obtained with 
the unitary transformation. The correct way to achieve the 
Hamiltonian He (p) is by the sub-space average of H(p) such 
that Hc{p) = {a\H(p) |a), where |a) is the optical coherent 
state with the usual definition a\a) = a\a). Then one is able 
to obtain the semiclassical ground state of Hc{p) with the 
standard variational method, in which the coherent state ja) 
acts as a trial wave function. Once doing so the BP emerges 
in both the Rabi and JC models independent of the RWA as it 
should be. 


B. Variational ground-state of the Rabi model and 
nonvanishing BP 

The proper semiclassical approximation begins with the av¬ 
erage of the original Rabi-mode Hamiltonian Hq (for Xnr = 
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A) in the coherent state |j^ |a), which leads to the following 
effective spin Hamiltonian 

He{a) = {a\Ho\a) 

= L0\a\^ +'^(7 z + \{a + a*)ax- (9) 

Its eigenvalues can be obtained by solving the energy eigen¬ 
value equation, 


H,{a)\^)=E{a)\i^). (10) 


The average energy 


E± (a) 


(V^±|Tfe(a) IV'i) 

uj \a\^ ± (a + a*)^ (11) 


is a function of a ,which is considered as a variational param¬ 
eter to be determined by the variation procedure, here \'tjj±) 
denotes the two eigenstates of the effective spin Hamiltonian 
He{a). Consequently the variational ground-state energy can 
be obtained by minimizing the energy function E± (a). It is 
easy to find that only E- (a) leads to the true ground state and 
the extremum condition is 


dE_{a) 

da 


( 12 ) 


which gives rise to the average photon number of the semi- 
classical ground-state 


0 , 


A < Ae 


v'AVw^-z/VlSA^, A>Ac 

and the average energy 

A<A, 


Egs — 


A>A, 

u) 16X^ ’ 


where 


( 13 ) 


(14) 


Ac = -V^ 


is the well known critical value of the quantum phase tran¬ 
sition from the normal to the superradiant phases in the N- 
atom Dicke model Si. Therefore the desired semiclassical 
ground-state reads 


lx) 


lags) \^-) 

I | 0 ) Iff), A < Ac 

1 lo^gs) (C- Iff) + |e)), A > Ac 


(15) 


with C± = =F a;z^/2/2A, and the relevant BP is calcu¬ 

lated as 


Igs 


0, A < Ac 
27r |ags|^ , A > Ac 


(16) 


Obviously the BP of the semiclassical ground-state in Rabi 
model is non-zero when A > Ac- For a JC model the re¬ 
sults are the same except a possible shift of the critical point 
Ac. The variational result can be verified in a more general 
sense by numerical diagonalization of the Hamiltonian iTp 
with Atva = A. The numerical values of BP as a function 
of the dimensionless coupling constant g = Xjuj are plotted 
in Fig.l. It may be worthwhile to remark that the average 
photon-number ags obtained from the variational treatment 
with the coherent state |a) may have a certain amount of de¬ 
viation with respect to the accurate results of full quantum 
mechanical formulation iSIl , and also to the numerical diag¬ 
onalization shown in Fig. 1. After all, the variational result in 
the coherent state is just a semiclassical approximation. 



FIG. 1: BP of the ground state in the Rabi model versus the dimen¬ 
sionless coupling parameter g = A/o; for different atomic eigenfre- 
quency v and field frequency lo. 


IV. GROUND-STATE BP OF A THREE-LEVEL ATOM IN 
THE QUANTIZED LIGHT-FIELD 


We now generalize the above BP formulation of two-level 
atom to a three-level atom in the quantized cavity field. With¬ 
out loss of the generality, we consider the A-type atom with 
the energy eigenvalues Ei = huji, for f = 1, 2, 3. The allowed 
dipole-transitions are assumed to be |3) o |2) and |3) O |1) 
but not for |1) o |2) in . Beyond the RWA, the Hamiltonian 
of this system reads 

= H^+r,{b+b^){\l) (3| + |3) (1| + |2) (3| + |3) (2|). (17) 

Here, 

E[q = ujQb^b -\- LOi |1) (1| -F tU2 |2) (2| -F |3) (3|, 

and ff is the coupling constant between atom and the field of 
frequency wq. b and b^ denote the corresponding field opera¬ 
tors. 

Similarly, the average of Hamiltonian E[^ in the field 
coherent-state |/3) (61/3) = (3 1/3)) becomes an effective 















4 


Hamiltonian of atomic operator only 

= (/3|i/^|/3) 

= ujq \j3^ + LOi \l) {l\ + 2mp(|1) (3| + 

I 

|3)(1| + |2)(3| + |3)(2|), (18) 

with the complex number f3 = u + iv considered as a vari¬ 
ational parameter to be determined. For simplicity, let us as¬ 
sume uji = lo 2 , and then the energy eigenvalue equation 

= (19) 

can be solved with the results given by 

E^{l3)=u;o\l3\^+oju 


and 

lfl|2 , “^1 + W3 , / (wi - ,022 

E±il3) =u;o\P\ + - - -±Y- - - + 8r]^uE 

One can easily check that the energy function E^{(i) is the 
lowest one i.e. E^{I3) < Eq{13), E^{I3). Therefore, the 
ground state energy can be determined from the extremum 
condition dE^{j5)/d(3 = 0 and the result is 




Wl, T] < F 

-2rf‘ uJo(uJi-uj3)^ I UJ1+UJ3 ^ Ti' 
~o -321^-^ 2 ’ " ^ ^ 

with the relevant variational parameter found as 

0, r] < F 

K - v>F 

327]^ ’ ' 


I 3 gs — 


( 20 ) 


( 21 ) 


where F = -^/wq (wi — 023 ) / 8 . Consequently, we apply a 
unitary transformation U{t) = exp(—iT6^6) to the Hamilto¬ 
nian F[^ and obtain 


H^(t) = U{t)H^U\t) 

= 7Jo^ + p(&e- + 6te—)(|1)(3| 

+ |3)(1| + |2)(3| + |3)(2|). (22) 


Thus, when r varies slowly from 0 to 27r, the semiclassical 
ground state = t/(r) I'i’gs) '''iH acquire a BP evalu¬ 

ated as 

Igs = * ^ dT (^) I ^ I i’gs (^) ) 

= 27r(V'g^| &+6|V’g,), (23) 

which leads to the result 


7^ 

Igs 


= 27r|/3g 


0 , 


T] < F 


77 >F 


3277^ 


(24) 


Again, the BP depends on the coupling constant rj as that in 
the two-level case. As a comparison the BP values are also 
evaluated by the numerical diagonalization of the Hamiltonian 
F[^. Fig. 2 shows the plots of ground-state BPs versus the 
dimensionless coupling constant g' — g/uiQ . It is seen that 
the numerical results are qualitatively in agreement with our 
semiclassical analysis. 



FIG. 2: Ground-state BP in the three-level atom system versus the di¬ 
mensionless atom-field coupling constant g' for oji = iJ 2 = —0.25, 


u)3 = 0.25. 


V. CONCLUSION AND DISCUSSION 

We present in this letter a general formulation of the BPs 
for both the JC model (with RWA) and the Rabi model (with¬ 
out the RWA). In the semiclassical approximation the varia¬ 
tional ground-state and the related BP for both Rabi and JC 
models are obtained analytically in a consistent manner. We 
argue that the vanishing BP of the Rabi model |[^ is due to 
the improper semiclassical approximation with simply replac¬ 
ing the bosonic operator by a complex C-number since the re¬ 
sulted Hamiltonian does not correspond to the original one in 
the semiclassical level. The valid effective spin Hamiltonian 
He{a) in the semiclassical approximation is achieved by the 
average in the field coherent state and thus the macroscopic 
(semiclassical) quantum state should be obtained by means of 
the standard variational method. We show in this letter that 
the BP for a generic Rabi model is indeed non-zero. This ob¬ 
servation is also been generalized to a three-level atom in the 
quantized cavity field. 

The BP in the variational ground-state, which is a displaced 
vacuum (i.e. coherent state), may possess a simple interpreta¬ 
tion as suggested by one Referee of the paper that the action of 
the unitary transformation Eq. (|2]) is to take this coherent state 
around the origin in phase space and thus the emerging BP is 
just the one of a harmonic oscillator. Following the Ref. 
the BPs beyond the RWA, obtained in the present work, should 
be also verified experimentally, in principle, with the usual 
cavity quantum electrodynamic system. The analogous ex¬ 
periment to measure the BP has been also designed in a single 
solid-state spin-quibt ll2ll] . 
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